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Abstract 

The local composite operator is analysed within the algebraic renormalization in 
Yang-Mills theories in linear covariant gauges. We establish that it is multiplicatively renor- 
malizable to all orders of perturbation theory. Its anomalous dimension is computed to 
two-loops in the MS scheme. 
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1 Introduction 



The possibihty that gluons might acquire a mass through a dynamical mass generation mech- 
anism is receiving increasing attention, both from the theoretical point of view as well as from 
lattice simulations. Effective gluon masses have been reported in a rather large number of gauges 
pp. For instance, the relevance of the local operator A^A""^ for Yang- Mills theory in the Landau 
gauge has been emphasized by several authors [21 El- That this operator has a special meaning 
in the Landau gauge follows by observing that, due to the transversality condition d^A""^ = 0, 
the integrated mass dimension two operator / d'^xA'^A'^^ is gauge invariant. Remarkably, the 
operator A'^A"'^ in the Landau gauge is multiplicatively renormalizable |^ E], its anomalous 
dimension being given El by a combination of the gauge beta function, /3(a), and of the 
anomalous dimension, 7^ (a), of the gauge field, according to the relation 

7A2(a) = - + 7A(a) I , " = ' 

Moreover, lattice simulations |3_ have provided strong indications of the existence of the conden- 
sate {^^^A"'^^, which is deeply related to the dynamical gluon mass. A renormalizable effective 
potential for this condensate in pure Yang-Mills theory has been constructed and evaluated in 
analytic form up to two-loop order in resulting in an effective gluon mass rrigiuon ~ 500Mey. 
The inclusion of massless quarks has been recently worked out 0. Another analytic study 
of {^A^yl"^^ can be found in 8 . Also, lattice simulations [01 of the gluon propagator in the 
Landau gauge have reported a gluon mass Wgiuon ~ 600Afel^. Concerning other gauges, an 
effective gluon mass has been reported in lattice simulations in the Laplacian ^01 and max- 
imal abelian |1H I12j gauges. It is worth underlining that the local operator A^A""^ of the 
Landau gauge can be generalized ^Sl to the maximal abelian gauge, which is a renormalizable 
gauge in the continuum El El • It turns out in fact that the integrated mixed gluon-ghost 
operator* / d^x (^^^^"/^ + Cc"c") is BRST invariant on-shell |13j . a property which ensures 
the multiplicative renormalizability to all orders of perturbation theory jl7| of the local 
operator (^^A'^A'^^ + ^c°c°^ . The analytic evaluation of the effective potential for the conden- 
sate + ^c"c"^ has not yet been worked out. Nevertheless, we expect a nonvanishing 
value for this condensate, which would result in a dynamical gluon mass. This is supported 
by the fact that a renormalizable effective potential for the mixed gluon-ghost operator has 
been obtained JHl in the nonlinear Curci-Ferrari gauge, yielding a nonvanishing condensate 
^^A^ + ^ccj, which provides a dynamical mass for the gluons. The Curci-Ferrari gauge shares 
a close similarity with the maximal abelian gauge. We expect thus that something similar should 
happen in this gauge. 

A gluon condensate ^A^yl"^^ has also been introduced in the Coulomb gauge 1201 order 
to obtain estimates for the glueball spectrum. Older works [22 (221 12H1 12^ already discussed 
the pairing of gluons in connection with a mass generation, as a result of the instability of the 
perturbative Yang-Mills vacuum. Also, the dynamical mass generation for the gluons is a part 
of the Kugo-Ojima criterion for color confinement |251 126j . See |2Z1 for a review. 

In this work we analyse the ultraviolet properties of the local composite operator A^^"^ in 
the linear covariant gauges, whose gauge fixing term is 

J d^x (b^df.A^^' + ^b%^ + ifd^^D^^c^ ^ , (2) 

*In the case of the maximal abeUan gauge the group index a labels the off-diagonal generators of SU{N), 
with a — 1, N{N — 1). The parameter ^ is the gauge fixing parameter of the maximal abelian gauge. 
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where 6" stands for the Lagrange multipher and a is the gauge parameter. Our aim is that 
of estabhshing some necessary requirements in order to study the possible condensation of this 
operator, which would imply the occurrence of dynamical mass generation in these gauges. 
Notice that, unlike the case of the Landau and maximal abelian gauges, the quantity / df^xA^^ 
is now not BRST invariant on-shell. However, we shall be able to prove that the local operator 
A^^ is multiplicatively renormalizable to all orders of perturbation theory. There is a simple 
understanding of this property. In linear covariant gauges, due to the additional shift symmetry 
of the antighost, i.e. c ^ c + const., the operator Aj^ doesn't mix with the other local dimension 
two composite ghost operator cc, which cannot show up due to the above symmetry. We remark 
that the renormalizability of Aj^ is the first step towards the construction of a renormalizable 
effective potential in order to study the possible condensation of this operator and the ensuing 
dynamical mass generation. 

The work is organized as follows. In Sect. 2 we derive the Ward identities for Yang-Mills 
theory in linear covariant gauges in the presence of the local operator Aj^. These identities 
turn out to ensure the multiplicative renormalizability of Aj^. In Sect. 3 the explicit two- loop 
calculation of the anomalous dimension of Aj^ is presented. 



2 Algebraic proof of the renormalizability of the local operator 

We begin by recalling the expression of the pure Yang-Mills action in the linear covariant gauges 
S = Sym + Sgf+fp (3) 



(fxF^^F''^'^ + I (Tx ( 6"9^.4"^ + -1,%" + lf^^'D^^^ 



where 

Df^d.S'^^-gr^'A^^. (4) 

In order to study the local composite operator A'^A"'^, we introduce it in the action by means 
of a BRST doublet of external sources (J, A), namely 

Sj = s j d^x QaA^A"^' + |aj) = j d^x Q JA^l^"^ + XA'^^d^d' + I , (5) 
where s denotes the BRST nilpotent operator acting as 



sa;, = 








sc" = 




slf = 


b\ 


sb" = 


0, 


sX = 


J, 


sJ = 


. 



(6) 

According to the local composite operators technique jSH]) the dimensionless parameter ^ is 
needed to account for the divergences present in the vacuum Green function (^A^(x)A^{y)), 
which turn out to be proportional to J^. As is apparent from expressions and Q, the action 
{Sym + Sgf+fp + Sj) is BRST invariant 

s {Sym + Sgf+fp + Sj) = . (7) 
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2.1 Ward Identities 



In order to translate the BRST invariance © into the corresponding Slavnov- Taylor identity 
[28j . we introduce two further external sources $7^, coupled to the non-linear BRST variations 



of A"', and 



with 



sQt = sL'^ = . 



Therefore, the complete action 



S = Sym + Sgf+fp + Sj + S, 



ext 



obeys the following identities 
• The Slavnov- Taylor identity 
5(S) = 



The linear gauge-fixing condition 



The antighost equation 



(8) 
(9) 

(10) 
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(12) 



(13) 



Notice also that the additional shift symmetry in the antighost present in the linear covariant 
gauges 

c ^ c + const. (14) 

is automatically encoded in the antighost equation p3() . Indeed, integrating expression (|13|1 on 
space-time yields 



d^x— = 



(15) 



which expresses in a functional form the shift symmetry ()14() . Equations ()13() , ()15() imply that the 
antighost field can enter only through the combination + d'^c"'), forbidding the appearance 
of the counterterm c"c". As a consequence, the local operator A^j^A""^ does not mix with c"c" in 
linear a-gauges. 

Let us also display, for further use, the quantum numbers of all fields and sources entering 
the action S 





A, 


c 


c 


b 


A 


J 


n 


L 


dim . 


1 





2 


2 


2 


2 


3 


4 


gh — number 





1 


-1 





-1 





-1 


-2 



(16) 
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2.2 Algebraic characterization of the general local invariant counterterm 



In order to characterize the most general local invariant counterterm which can be freely added 
to all orders of perturbation theory [22] , we perturb the classical action S by adding an arbitrary 
integrated local polynomial J]^"""^* in the fields and external sources of dimension bounded by 
four and with zero ghost number, and we require that the perturbed action (S+eS™""^*) satisfies 
the same Ward identities and constraints as S to the first order in the perturbation parameter 
e, i.e. 



+ eS' 



count 



Ua 



df'A''^ + + 0(e 



2\ 



This amounts to impose the following conditions on J]"^"*^"^* 
and 

rvcount 

— + *'^ = °' <^°> 

where Bs is the nilpotent linearized operator 

/■ ,4 / (5S 5 6J: 6 (5S 5 6^ 6 ,^6 

(21) 

Ss0s = . (22) 

Taking into account that (J, A) form a BRST doublet, from the general results on the cohomology 
of Yang-Mills theories [311] it turns out that the external sources (J, A) can contribute only 
through terms which can be expressed as pure Ss- variations. It follows thus that the invariant 
local counterterm S'^""'^* can be parametrized as 



Tcount 



I d^xF^.F'^^^ + SsA-i , (23) 



where o" is a free parameter and A ^ is the most general local polynomial with dimension 4 and 
ghost number —1, given by 

(24) 

with ai, ,08 arbitrary parameters. From the conditions (|19() . H20|l it follows that 

0-3 = fli; 0.4 = = aj = , (25) 
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so that A ^ reduces to 

= J (fx (^ai{% + d^lf)A'''' + aiL^c" + ae^^'^^A^l + ag^AJ^ . (26) 

Notice that the vanishing of the coefficient 07 imphes the absence of the counterterm Jc"c". As 
aheady underhned, this ensures that the operator A"'^A'^ does not mix with the ghost operator 
c"c". Therefore, for the final form of the invariant counterterm one obtains: 

+^(2ai + a6)JA\A'"' + (ai + ae - a2)X^^,c''A''^' 

-^gfatcL'^c'c^ + ^^j') . (27) 

It remains to discuss the stability of the classical action |2Sj, i.e. to check that 1]'^°'^°* can be 
reabsorbed in the classical action T, by means of a multiplicative renormalization of the coupling 
constant g, the parameters a and ^, the fields {(j) = A, c,c, b} and the sources {$ = J, A, L, il}, 
namely 

S(5, e, a, (/>, ^) + = ^{go,Co, «o, (^0, ^0) + O(e') , (28) 

with the bare fields and parameters defined as 



a 


— 7^1'^ A°- 

— ^ii 


Qa _ 
) "OAt ~ 




, 50 


= Zgg , 


Co 


= zli'^e 


ra _ 
) ^0 — 


ZlL'' 


, ao 


= Zo,a 


Co 


= z]i^i^ 


, Jo = 


ZjJ 


, eo 


= z^i, 


IP' 


= zl'h'^ , 


Ao = 


ZxX . 







(29) 

The parameters a, ai, 02, ag, turn out to be related to the renormalization of the gauge 
coupling constant g, of c"^, J, A, and ^ , according to 

Zg = , 



ZT = l + e[^- + ai 

^1/2 ^ ^_^/^l+a2 
V 2 

Zj = 1 + e (ae - 0-) , 

7 1 , ( ai-a2-a 
Z\ = l+eloeH 

Z^ = 1 + e (as - 2a6 + 2a) . 

(30) 

Concerning the other fields and the sources 0^, L"", it can be verified that they are renormalized 
as 

Zc = Zc , Ziy = z^^ , Zq = zy^ 

Zl = Z'l^ , Z^ = Za. (31) 



2 



cr 
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This completes the proof of the multiplicative renormalizability of the local composite operator 
Aj^ in linear covariant gauges. Finally, it is useful to observe that, from eqs.ljSOJ, one has 

= Zj Zy^z'/, (32) 

from which it follows that the anomalous dimension of A^, turns out to be related to that of the 
composite operator A^j^d^d^ 

lAdc = lA^ +Ic + 1A , (33) 

where 7c, 7^, 7^2, and 7a<9c are the anomalous dimensions of the Faddeev-Popov ghost c", of 
the gauge field j4^, of the operator A^, and of the composite operator A'^d'^c'^, which are defined 

as 

7c = i^d^ In Z]/^ , 7A = fid^ In Z^J^ , 7^2 = fid^. In Zj , -fAdc = l^d^ In Z^ , 

(34) 

where /i is the renormahzation scale. As expected, property ()33() relies on the fact that A'^d^c"' 
is the BRST variation of i-e. 

s^-^ = -Ald^^c'^ . (35) 

Although we did not consider matter fields in the previous analysis, it can be checked that the 
renormalizability of and the relation H33|l remain unchanged if matter fields are included. 



3 Calculation of the two-loop anomalous dimension of A^^ 

We now turn to the computation of the anomalous dimension of in an arbitrary linear gauge. 
The method exploits the lack of mixing in the linear covariant gauges between A^^ and the other 
dimension two Lorentz scalar zero ghost number operator c"c°, which we have already noted. 
For instance, in the Curci-Ferrari gauge although both operators mix there is a combination, 
O = \A?^ + acf'c"', which remains multiphcatively renormalizable. Prior to the proof of J, that 
the anomalous dimension of O was related to the /3-function and the gluon anomalous dimension, 
70 (a) was explicitly computed at three loops in MS in |S] . That method involved substituting the 
operator in a ghost two-point function with a non-zero momentum flowing through the operator 
itself and one external ghost momentum nullified. This configuration allowed for the application 
of the Mincer algorithm, written in the symbolic manipulation language Form, |32[ I33j . 
A ghost two-point function was chosen to avoid the appearance of spurious infrared infinities 
which would arise for this momentum configuration if the external legs were gluons. To determine 
7^2 (a) in the linear gauges we are forced into the same approach as [Sj due to the infrared issue 
with gluon external legs. Hence, we have renormalized the momentum space Green's function 
{c''{p)\[Al]{-p)c^{ld)) where p is the external momentum. Clearly, this has no tree term and 
therefore to deduce 7^2(0) at n-loops requires renormalizing the Green's function at (n+l)-loops 
as the one loop term corresponds to the tree term of (^^(p)^[^^](— p)^^(O)). This is evident, 
for example, by drawing one and two loop diagrams for the various Green's functions based on 
the interactions of the Yang-Mills action, eq.(|21). Since the Mincer algorithm currently only 
extracts the simple poles in e in dimensional regularization to three loops, where d = 4 — 2e, 
this means we have only computed 7^^2(0) to two loops. Though this will be sufficient to deduce 
the effective potential of A^^ to one loop. The Feynman diagrams for our Green's function are 
generated with Qgraf, |34j . and converted into Form input notation, At one loop there is 
one diagram which plays the role of the tree diagram and at two loops there are 15 diagrams. 
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The bulk of the calculation, however, is in the evaluation of the 314 three loop graphs. Since 
there is no operator mixing we can apply the rescaling technique of jSS] for automatic multiloop 
computations to find the renormalization constant Zj^2 . From this we deduce 



7a2 (a) 



[(35 + 3a) Ca - l^TpNf] 
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(449 



+ 33a + ISa^ ) C\ - ISOCATpNf 



192CFTFNf\ — + 0{a^) (36) 



in the MS scheme where is the number of quarks and the colour group Casimirs are defined by 
Tr (t'^T*) = T°-T'^ = Cpl and /«crfj6cd ^ c^^afe. in deriving ^ from the corresponding 

renormalization constant we have verified that the double pole in e is correctly reproduced for 
all a. Moreover, reduces to the Landau gauge expression of [H^ . 



4 Conclusion 

We have investigated the renormalizability of the dimension two operator Aj^ in arbitrary co- 
variant linear gauges in Yang- Mills theories, due to the possibility that it might condense and 
develop a non-zero vacuum expectation value. This would generalize to these gauges previous 
results obtained in the Landau gauge |2l 13 IH [3 . One feature of our analysis is that, unlike the 
Curci- Ferrari gauges the operator does not mix with the other dimension two local 

composite operator c"c°. This is a general feature of the linear covariant a-gauges, present also 
in the Landau gauge jBj, a = 0, which is a consequence of the additional shift symmetry in 
the antighost (|14|) . Importantly the operator Aj^ can thus be treated in isolation as it does not 
require any ghost dependent operator. 

Finally, we underline that the multiplicative renormalizability of Aj^ and the explicit knowl- 
edge of its anomalous dimension for all a, eq. (|36|l . are central ingredients towards the con- 
struction of a renormalizable effective potential for studying the possible condensation of this 
operator and the related dynamical mass generation, as was carried out in the Landau jl||7j and 
Curci-Ferrari gauges. 
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